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Abstract 

A siraple molecular thermodynamic approach is applied to the study of the 
adsorption of gases of chain molecules on solid surfaces. We use a model based 
on the Statistical Associating Fluid Theory for Variable Range (SAFT-VR) 
potentials [A. Gil-Villegas, A. Galindo, P. J. Whitehead, S. J. Mills, G. Jackson, 
A. N. Burgess, J. Chem. Phys. 106 (1997) 4168] that we extend by including 
a quasi-two-dimensional approximation to describe the adsorption properties of 
this type of real gases [A. Martinez, M. Castro, C. McCabe, A. Gil-Villegas, J. 
Chem. Phys. 126 (2007) 074707]. The model is applied to ethane, ethylene, 
propane, and carbon dioxide adsorbed on activated carbon and silica gel, which 
are porous media of significant industrial interest. We show that the adsorption 
isotherms obtained by means of the present SAFT-VR modeling are in fair 
agreement with the experimental results provided in the literature. 
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1. Introduction 



The main component of natural gas is methane, but small amounts of other 
gases such as ethane, ethylene, propane, propadiene, butane, methylpropane, 
2-methylpropane, and carbon dioxide are present. These components are in 
negligible amounts in natural gas but become a problem when the methane is 
stored through adsorption and are less volatile than the methane and therefore 
preferentially occupy the active sites of the adsorbent. To solve the problem 
of this type of pollutants one must know their specific adsorbate properties on 
various adsorbents at ambient temperatures and a broad range of pressures. 

The study of adsorption is very important for designing a good adsorbent. 
A correct interpretation and quantification of the adsorption isotherms are re- 
quired for the development of adsorption technologies in the case of the separa- 
tion and purification of gases and liquids [H 01 ■ A detailed molecular thermo- 
dynamic model based on perturbation theory for fiuids with highly anisotropic 
interactions has been developed by Wertheim [5| some time ago and has been 
known as the Statistical Associating Fluid Theory (SAFT) since the works of 
Chapman and collaborators 0, [sj . 

The SAFT approach in three dimensions has been used in a variety of sys- 
tems and complexes relevant to the industry, and due to its predictive power 
and versatility started to be widely used in chemical engineering. Recently, 
Martinez and collaborators [ij modified the two dimensional form of SAFT-VR 
approach (SAFT-VR-2D) to a quasi two-dimensional variant and obtained very 
good results in the case of adsorbed gases such as methane and nitrogen on ac- 
tivated carbon and in addition for as complex molecules as asphaltene adsorbed 
on some natural porous rocks present in an oil well 

This paper presents the application of the SAFT-VR formalism, both in the 
context of 2D p, 0, S] and 3D ^ , to obtain the adsorption isotherms of ethane 
{C2H4), ethylene {C2HQ), propane {C^Hs), and carbon dioxide (CO2) in the 
case of activated carbon and silica gel as adsorbents at pressures up to 4 MPa 
and for three different temperatures. 

2. Theory 

In this section, we summarize the principal terms that enter the SAFT-VR 
model in 2D and 3D and the way we combine them to get the quasi two- 
dimensional approximation. 

2.1. Adsorption of monomeric fluids 

The model consists of a sirnple fiuid of N particles with diameter a in the 
presence of a uniform wall 0] . The interaction exerted by the wall on the 
particles is of the square well (SW) type 




z < 

< z < Xwcr 
X^a < z , 



(1) 
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where z is the perpendicular distance of the particles from the wall, e^u is the 
depth of the well, and \w<J is the range of the attractive potential. 

The system is divided into two subsystems: a fluid whose molecules are near 
the wall forming the adsorbed fluid considered as a substrate of the wall and a 
fluid whose molecules are far from the wall forming the fluid in the bulk. The 
fluids, both the adsorbed and the bulk phases, have different properties because 
of the presence of the wall. In particular, it is well known that the interaction 
between the molecules depends on how close they are to the wall, in other words 
the binary interaction between particles is different in the adsorbed and the bulk 
phases Thus, the particle density in the fluid, p, is a function of z. To 

describe the amount of adsorbed particles, we define the following parameter 

/•oo 

r = / \p{z) - pk]dz , (2) 



where pi, is the particle density in the bulk region. Since the length scale in the 
adsorbed fluid is defined by Au,cr, we can rewrite Eq. ([2]) as 

r = / p{z)dz - pbX^cr (3) 
Jo 

where the integral in the right hand side of Eq. ([3]) is the density of the adsorbed 
particles, Pads- In thermodynamic equilibrium, the chemical potential of the 
adsorbed phase Pads and the bulk phase /if, should be equal. Therefore, Pads 
can be obtained from the condition 

Pads = Pb ■ (4) 
The partition function of the adsorbed fiuid is given by 

K^ds 

Zads = ^|^37v Qads , (5) 

where Vj^^ is the volume containing the adsorbed fluid, A is the de Broglie ther- 
mal wavelength in terms of Planck's constant, and Qads is the configurational 
partition function of the adsorbed system. 

Introducing the adsorption area S, we have Vads = SXwCr, and Qads can be 
rewritten as 0, 0: @1 

Qads = Q1DQ2D , (6) 

where Qid and Q2D denote the configurational partition functions in one and 
two dimensions, respectively, and are given by 

Q2D^^J dx'^dy'' exp (^-p^i-l^{x,y)^ (7) 

and 

Qic -exp(-7V/3up^(z*)) , (8) 
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where z* is that value of the coordinate z that guarantees the mean value of the 
Boltzmann factor [gJ . 

Then the final expression for the partition function has the following form 

^AT = ^^z^^'Qfz, (^^) exp(iV/3upJ . (9) 

Applying the standard relation A — —kThiZ, the Helmholtz free energy of the 
adsorbed fluid is given by 

^ ^ _ In pu,^ , (10) 

NkT NkT V As / 

where A2D is the Helmholtz free energy of a two-dimensional fluid with interac- 
tions given by the binary potential 0(a;, y). In perturbation theory up to second 
order, the dimensionless Helmholtz free energy per particle, A2D / NkT , can be 
written as follows 

- ln(p<,,,A2) + + P'ar • (11) 



NkT ' NkT 

Here Ahd is the Helmholtz free energy for hard disks, (3 = 1/fcT, af^ and 
are the first two terms of the 2D perturbation expansion. It is worth mentioning 
that the expression for Qio as appears in Eq. ([5]) is accurate and can be used 
to define the energy parameter of the wall, e^,. For the simple case where the 
wall-particle interaction is given by a square well with range X^a and energy 
depth ew, we get Upw{z*) — — e^u. 

A similar expression in perturbations is considered for the fluid in the bulk 



Hi 



, , .3-. , 2D , n2 2D /nx 

= InipbA ) + + /3ai + (3 . (12) 

Getting the chemical potentials fiads and /i;, from Eqs. ([T0| . (fTTj) . and ([T2|) . we 
can rewrite the condition as follows 

Man = ^J'2D + IJ-io , (13) 

where and /12D are the chemical potentials in ZD and 2D, respectively, and 
is the wall contribution to the chemical potential 

Mu. = ln(A^,) + Up^{z*) . (14) 

2.2. Adsorption of chain molecules 

In the S AFT- VR framework the theory presented in Subsection 12.11 can be 
easily extended to the more general case of adsorption of homonuclear and het- 
eronuclear chain molecules [111] , which is one of the tasks of this work. In the 
simpler case of homonuclear molecules, we consider the fluid as made up of N 
molecules each composed of m identical spherical units (monomers) of diameter 
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(7. We assume that both the particle-particle and particle-wall interactions are 
described by square-well potentials. The thermodynamic equilibrium between 
the bulk phase and the adsorbed phase of the fluid is obtained by imposing the 
condition in Eq. (1131) of equal chemical potentials. The latter quantities are 
derived from the SAFT-VR expressions for the Helmholtz free energies of chain 
molecules in 3D and 2D. 



SAFT-VR approach. The Helmholtz free energy A of chain molecular gases 
is composed of three separate terms 0, : one refers to the ideal gas free energy, 
another to the free energy of monomeric units and finally the term corresponding 
to the formation of chain molecules incorporating a given number of monomeric 
units 

A j^ideal j^mono j^chain 

NkT ~ NkT ^ NkT ^ NkT ' ^^^^ 
We will now provide the main expressions of each term in this equation in 3D 
and 2D. For more details the reader is directed to refs. [1, 0, @]- 
The ideal contribution is given by: 

Aideal 
Aideal 

- MPa^sA') , (17) 

where A is the de Broglie thermal wavelength. The monomer contribution is 
given by: 

^ (18) 

NkT ^NskT ' 

A mono A 

^ (19) 

NkT NskT ' 

where N^ is the number of monomer units of the chain molecule, and A^d and 
A2D are the free energy of monomer fluids in 3D and 2D, respectively, that can 
be obtained from perturbation theory according to Eqs. (ITU)) and (fT^ . 
Finally, the chain contribution is given by: 

A chain 

^ = -(™-l)lny3^'^^(a), (20) 

A chain 



(™-l)lny^^^^^(a) , (21) 

where y^^'^^ and y^-D.SW g^j-g ^-j^g correlation functions in the 3D and 2D back- 
grounds obtained from the corresponding radial distribution functions through 
the well known relationship y{r) = g{r)e^"^'^\ Following the SAFT-VR proce- 
dure, the functions g{r) are obtained as perturbation expansions 

g^^^^^ia)^9"'{a)+f3eg!^{a), (22) 

g2AW(a)^.g^^(a)+/3eg?^(a), (23) 

where g^^ and g^^ are the radial distribution functions of the hard spheres 
and hard disks, respectively. 
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Table 1: The values of molecular parameters used to describe the adsorption of ethane, ethy- 
lene, propane and CO2 on activated carbon. 



Substance 



TO a{A) A e/fc(K) 



eadslk{K) 



/e T(K) 



C2H4 
C2HQ 

CO2 



1.2 3.743 1.444 236.1 

1.33 4.233 1.449 224.8 

1.67 3.8899 1.4537 260.91 

2.0 2.7864 1.5257 179.27 



1.1592 
1.1592 
1.16296 
1.2620 



188.88 
179.84 
208.728 
143.416 



0.2453 7.80 301.4 

0.2453 9.8 301.4 

0.8165 10.0 338.7 

0.8165 7.7 301.4 



3. Results 

A comparison of the theoretical adsorption isotherms with the experimen- 
tal results for the following gases C2H4, C2HG, C^Hs, and CO2 onto activated 
carbon and silica gel is displayed in Figs. [TlfT2l at different temperatures and 
pressures. Table I and II show the parameters we have used to describe the 
investigated simple fluids both in the bulk phase and the adsorbed phase within 



ethylene, and propane are represented as fluids of non-spherical monomers 
(to — 1.2,1.33,1.67), while CO2 is represented as a dimer fluid (to — 2). The 
molecular diameters a are the same both in the bulk and the adsorbed phases, 
while the parameters of the SW attractive potential are different from case to 
case. These differences are due to the influence of the wall on the substrate of 
particles. As mentioned in previous works @, 0, @Ij according to the theoreti- 
cal results derived by Sinanoglu and Pitzer for a Lennard- Jones fluid 10], the 
energy of the potential well depth for the particle-particle interaction in the ad- 
sorbed monolayer is reduced by a factor of 20%-40% with respect to their bulk 
phase values. We chose a reduction of 20% in all cases. Once eads is selected 
in this way, the range Xads can be determined by reproducing the experimental 
ratios between the critical temperatures of the bulk phase, Tj^, and the adsorbed 
phase, T"'*'', i.e., Rc = T^"^^ /T^. The value of this ratio is known for the case 
of noble gases and methane adsorbed on graphite surfaces, Rc ~ 0.4. Next, 
Aads was determined from this ratio and the values of the set of parameters 
TO, a, e. A, tads for ethane, ethylene, propane and carbon dioxide are those re- 
ported in Table 1 and 2. We notice here that at higher confinement regimes 
for the fluid, Rc could depend on the degree of confinement for some types of 
activated carbon as well as for silica gels and therefore, strictly speaking, our 
results apply in the low confinement regime. 

For the range parameter A^, we used two values, 0.2453 and 0.8165. The first 
one corresponds to the value found for the adsorption of noble gases on graphite, 
whereas the latter one is the upper limit used to describe the adsorption of a 
monolayer according to the mean- field approach [l6| . The energy parameter was 
adjusted to reproduce the experimental adsorption isotherms. In Figures [TJH 
the reported results for the adsorption on activated carbon were conducted at 
temperatures of 301.4 and 338.7 K at pressures up to 3.5 MPa. One can notice 



the SAFT-VR approximation 




6 



Table 2: The values of molecular parameters used to describe the adsorption of ethane, ethy- 
lene, propane and CO2 on silica gel. 



Substance 


m 


a{A) 


A 


e/fc(K) 


Aads 


eads/k{K) 


Au, 




T(K) 


C2H4 


1.2 


3.743 


1.444 


236.1 


1.1592 


188.88 


0.2453 


4.5 


278.0 


C2H4 


1.2 


3.743 


1.444 


236.1 


1.1592 


188.88 


0.2453 


3.5 


303.0 


C2HQ 


1.33 


4.233 


1.449 


224.8 


1.1592 


179.84 


0.2453 


6.3 


278.0 


C2HQ 


1.33 


4.233 


1.449 


224.8 


1.1592 


179.84 


0.2453 


5.0 


303.0 


CsHs 


1.67 


3.8899 


1.4537 


260.91 


1.16296 


208.728 


0.8165 


2.69 


293.0 


C3HS 


1.67 


3.8899 


1.4537 


260.91 


1.16296 


208.728 


0.8165 


3.0 


303.0 


CO2 


2.0 


2.7864 


1.5257 


179.27 


1.2620 


143.416 


0.8165 


4.95 


293.0 


CO2 


2.0 


2.7864 


1.5257 


179.27 


1.2620 


143.416 


0.8165 


6.0 


303.0 



the good agreement between theory and experimental data, unless in Figure [3] 
where our SAFT-VR calculations although not very close to the experimental 
data still follow their main trend. Figures [Mill show the results of adsorption 
on silica gel, for three different temperatures, i.e., 278, 293, and 303 K. We can 
see from these figures that the theoretical behavior agrees with the experimen- 
tal data, except at high pressures where a deviation is noted which may reflect 
the fact that the molecules are not spherical. Moreover, in the case of carbon 
dioxide the quadrupolar contribution is not taken into account and it is known 
that it influences the thermodynamic properties of both the bulk phase and the 
adsorbed phases. Another feature that is not taken into account in the model is 
the way the molecules adhere to the wall, that is, if their chain structure is per- 
pendicular, parallel or diagonal to the wall. Despite these limitations, the model 
provides very good theoretical results compared to the available experimental 
data. 

An advantage of this SAFT-VR thermodynamic model is that one only needs 
to know the particle-particle and the wall-particle interaction energies together 
with the molecular diameter, and requires the adjustment of only one parameter 
that corresponds to the strength of the attraction energy of the wall, but this 
adjustment is based on the experimental data of the adsorbate energy 11^,, which 
is known in the literature as the isosteric heat of adsorption. For example, in the 
case of krypton adsorbed on graphite, C/^=ll-72 kJ/mol [2^ and corresponds 
to a value of e=8.168 [16|, which is of comparable magnitude to the values used 
in this work. On the contrary, other models use at least four parameters to 
make adjustments and to reproduce the experimental data of the adsorption 
isotherms and usually these models are semi-empirical 21 1. 



4. Conclusions 

We have applied the SAFT-VR approach to model adsorption isotherms 
of simple fluids with molecules of non spherical shape. To achieve this goal 
we used a SAFT-VR approach in which the chemical equilibrium condition in 
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both 3D and 2D has been employed. For the cases of real fluids considered it 
is demonstrated that a good fit of these theoretical predictions is obtained by 
adjusting only a single parameter, which is the energy of the attraction interac- 
tion of the molecules towards the wall. We improved further the predictions by 
taking into account the positions of the particles on the wall as well as adding 
a new term to the SAFT-VR-2D corresponding to the quasi-two-dimensional 
quadrupole approximation. We recall that other semiempirical models available 
in the literature [2^ make use of at least four adjustable parameters with- 
out a clear physical interpretation, whereas the SAFT-VR parameters have a 
molecular basis being related to the molecular forces acting in the adsorption 
process. 

On the other hand, the complicated porous structure of silica gels and some 
of the activated carbons is at the origin of the fractal character of their surfaces. 
Taking into account this fractality [^l and the distribution of pores in the 
SAFT-VR approach implies introducing fractal type parameters that remains 
as an interesting issue for future investigation. 
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Figure 3: Adsorption of CsHg on activated carbon at 338.7 K. Tlie circles are experimental 
data from ref. [18l |. 
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Figure 4: Same as in Fig.[3l but for CO2. 
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Figure 5: Adsorption of C2H4 on silica gel at 278 K. The circles are experimental data from 
ref. [H 
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